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Abstract 

s 

Clean metallic superlattice systems composed of alternating layers of super- 

o ' 

O ■ conducting and normal materials are considered, particularly aspects of the 



X 



proximity effect as it affects the critical temperature. A simple model is used 



to address the question of when a finite-sized system theoretically approx- 
imates well a true infinite superlattice. The methods used in the analysis 
afford some tests of the approximation used that the pair amplitude of the 
Cooper pairs is constant over a superconducting region. We also use these 
methods to construct a model of a single superconducting layer which in- 
tends to incorporate a more realistic form of the pair amplitude than a simple 
constant. 
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I. INTRODUCTION 



Metallic superlattice systems consisting of alternating layers of superconducting and 
normal films have attracted wide interest, both experimentally and theoretically [|IJ. Much 
of the theoretical work has gone into the understanding of the proximity effect in such 
systems. This effect between adjacent layers of materials with different superconducting 
properties has been studied extensively for the case of two such neighbouring films. Close 
to the transition point Ginzburg-Landau theory can be employed [2-7], but away from 
this point other approaches such as the constant pair amplitude approximation [8-15], a 
"Cooper" proximity effect argument [ l6| , the tunneling model of McMillan [17], a WKB 
approximation or a semi-classical approximation [l9[j can be used. 

Implicit in these approaches when applied to superlattices is the argument that, from 
symmetry considerations, the form of the pair amplitude of the Cooper pairs for a bilayer 
is the same as that for a superlattice consisting of an infinite number of adjoining bilayers. 
While true for an infinite superlattice, experimentally these systems are of a finite extent, 
and one could then question the reasonableness of the bilayer assumption for finite systems. 
This issue could also arise in the general context of the applicability of Bloch's theorem for 
periodic potentials in cases of large but finite-sized systems, as for instance as applied to the 
Kronig-Penney model |2"0fl . In this paper we will address some aspects of this question of 
how many bilayers are effectively needed to have an infinite superlattice system, particularly 
as related to how the critical temperature of such systems is affected. 

The paper is organized as follows. To introduce the approach and to establish notation, 
in Sec. II we review a method of analyzing the proximity effect based on the constant 
pair amplitude approximation as applied to the simple case of a single superconducting film 
immersed in a normal metal [|U|]. In Sec. Ill we extend this method to an infinite superlattice 
composed of alternating superconducting and normal layers. To test some of the ideas we 
shall later employ, we consider the superlattice as the limit of a finite sized system as the 
number of layers approaches infinity; the result we find for the critical temperature agrees 
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with that obtained by arguments which explicitly incorporate Bloch's theorem using the 
same general approach [^TJ. In Sec. IV we consider finite-sized "superlattices" . We argue 
from the calculated form of the pair amplitude that the approximation of a constant pair 
potential is reasonable in two distinct cases: for thick superconducting films over a wide 
range of thicknesses of the normal film, and for thin superconducting and normal films. 
As related to finding the critical temperature in this approach, in the first case a relatively 
small number of layers are needed before the finite-sized system approximates well an infinite 
superlattice, whereas in the second case quantitative differences arise between the finite and 
infinite lattices. In Sec. V we reconsider the single superconducting layer of Sec. II regarding 
the issue raised in the previous section of the reasonableness of the constant pair potential 
assumption. We study a model based on these methods which aims to incorporate in a 
simple manner a more realistic form of the pair potential, and find that the resulting zero 
temperature critical thickness gets enhanced compared to that of Sec. II found using the 
constant pair potential approximation. Sec. VI contains some concluding remarks. 



II. A SINGLE FILM 



We first study the case of a single superconducting film embedded in a normal metal [|T5 
The film, of width 2a, is centered at the origin x = 0, as in Fig. 1. The Hamiltonian for this 
system in terms of the electron wave function ip is given in the mean field approximation by 
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H = I d 3 rg{x) | F(x) | 2 + J rf 3 r0 f (r) 
- / d 3 r g(x)F(x)(^(r)T 1 (j)(r), 
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Ts0(r) 
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where Tj are the Pauli matrices, F(x) = F*(x) = — (^(x^^x)) is the pair amplitude, 




(2) 



and g(x) denotes the spatially inhomogeneous BCS coupling constant. We assume the 
superconducting and normal metals are clean and have the same effective electron masses, 



Fermi velocities, Debye temperatures, etc. 
The Nambu doublet is then decomposed as 

1 



*) = Et^i/ dE J dH [u {i) (E,l,x)a^(E,l)e-^ m - hp) 



(27T) 

+ v {i \E,l,x)b {i) \E,l)e i{Et ~ vp) ] 



(3) 



where 1 = (k y ,k z ), p = (y,z), and "i" labels the particular linearly independent solution 
present. If we now introduce 



u 



V(E,l,x)= X {t \E,l,x)e^, 



(4) 



where q = ykp — l 2 , then assuming x(E,l,x) to be smooth on the atomic scale results in 
the Bogoliubov equation following from Eq. ([[]): 



E + i^i- + g(x)F(x)r 1 ) X (i) (E, I, x) = 0. 



(5) 



The corresponding wave function v(E, I, x) = u(—E, I, x) can then be readily obtained. 

The non-linearity introduced by F(x) = —{ip^(x)ipi(x)) makes Eq.(H) difficult to solve for 
inhomogeneous geometries. To obtain a simple solution we therefore make the approximation 

A for —a < x < a 

(6) 

for x < —a or x > a 
where A is a constant. The solutions to Eq.(|5|) under this approximation are readily found 
to be 
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1 + 7 



for x < —a or x > a, 



e ipx for —a < x < a 



(7) 



where 



p = jk = mEj/q = mEj/ \jk 2 F — l 2 , 



5 = -A/E, 
7 = v 7 ! 3 ^ 2 . 



(8) 



The coefficients of the wave function entering into the solutions of Eq. (|7|) can be determined 
by imposing continuity at the interface x = ±a together with the normalization conditions 

/OO pCsO 
dxu(E, I, x)u{E', I, x) = 2tc5(E — E') = / dxv(E,l,x)v{E' ,l,x), 
-oo J —oo 

dxu(E,l,x)v(E',l,x) = 0. (9) 



Thus far these solutions contain the free parameter A. We fix this parameter by requiring 
that it minimize the free energy of the system: 

dH 

where the angular brackets denote the thermal average. Using Eq.(|l]) we find that this 
condition translates into the self-consistent gap equation 

A= |- r F(x,A)dx. (11) 

ZCL J —a 

In other words, we demand that the spatial average of the calculated pair amplitude be 
equal to the constant A assumed in finding the solution. 

In the limit A — > the gap equation ( |TTj ) will determine the critical temperature T c (a) as 
a function of the thickness a of the superconducting film. Using the solutions and boundary 
conditions discussed previously, we find this equation becomes, assuming T c (a) <C T d , 

I R roo 

— - = H{u d ) - | j Q duH(u) cosh- 2 (^/2), (12) 

where ft = 1/ksT, N(0) = mk F /2ir 2 is the density of states of one spin projection at the 
Fermi surface and, with A = and Q = u/Ep, the function H(u) = H(AQ) is given by 

H(x) = 7 e + ln(x) - -G(x), 

G(x) = cos(x) - + 2 Ci(x) + x si(x). (13) 



Here Ci(x) and si(x) are the Cosine and Sine integrals defined as p2 

r x sin t 



Si(x) = si(x) + — = / 
2 Jo 



-dt, 



t 

rx QQg £ — \ 

Ci(x) = 7 e + In x + / dt, (14) 
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and 7 e is Euler's constant. The function H(x) of Eq.(|T3|) has the following limiting forms: 



lim H(x) = 7 e + ln(x) + 0(1/ x) 



TTX 



lim H(x) 



0{x 2 ). 



(15) 



For comparison with other geometries the numerical solution to Eq.(|12|) giving the critical 
temperature T c (a) will be discussed in Sec. IV. For now we only note that Eq.(|l!J) can be 
solved analytically in the zero temperature limit, which gives the critical thickness a c below 
which superconductivity cannot be maintained at any temperature: T c (a c ) = 0. One finds, 
assuming A c Q d 3> 1, 

1 fTp\ 1 



-7re~ 7e £ ~ 0.882^ 



(16) 



7rkp \T C J 2 

where T c is the bulk critical temperature and £o = vf/tcAq is the bulk coherence length. This 
value of the critical thickness compares favourably with that found using other approaches 
|7|JT7[|, and will be used in subsequent sections as a convenient length scale when discussing 
other geometries. It is also useful for later comparisons with other geometries to introduce 
a c directly into the gap equation of Eq.(|T2"D; doing so leads to the following equivalent forms: 
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(18) 



The feature of these and the following gap equations that the thicknesses of the supercon- 
ducting and normal layers scale with a c makes it possible that dirty materials might also be 
described within the present framework by a simple reinterpretation of a c , as was found in 
the single layer case by a comparison with the tunneling model [|I"5l ,|l7|l . 
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III. INFINITE SUPERLATTICE 



We next study an infinite superlattice consisting of alternating superconducting films 
of width 2a separated by normal films of width 26. This is illustrated in Fig. 2. Using 
the approach of the previous section, we first consider iV such films and eventually take 
iV — > oo. We assume the pair potential g(x)F(x) in each superconducting layer is a constant 
A and vanishes in each normal layer. This allows us to solve the corresponding Bogoliubov 
equations (|5|) in each region subject to the appropriate boundary conditions. The parameter 
A is then determined by the minimization of the free energy of the system. This condition, 
in the limit A — > 0, gives the critical temperature T c (a, 6) as a function of the thicknesses a 
and 6 of the films. Taking the limit N —>■ oo then results in the gap equation 

\n(a/a c ) = I(r) - J(r,y), (19) 



where y = 1 + b/a, r and I(x) appear in Eq. (|18l) , and 

f°° dtsm 2 (t) , . , , . . . 
J(x,y) = xJ ^A^cot[ty]lncosh(t/x). (20) 

For infinite b the function J(r,y) vanishes, resulting in the gap equation of Eq.(^) for 
the single layer. J(r, y) thus represents coherence effects the other layers introduce to the 
single film geometry. Eq. (|T9"D agrees with the results of an analysis directly incorporating 



Bloch's theorem in this context for periodic systems | 21| . We present in the next section 
the numerical solution to Eq . (|19D in various cases and compare these results with those 
corresponding to a finite number of layers using this same general approach. 



IV. FINITE SUPERLATTICE 

In this section we examine by these methods a "superlattice" with a finite number of 
layers. As with the true infinite superlattice this consists of alternating superconducting 
films of width 2a separated by normal films of width 26, as in Fig. 2, but now after the last 
superconducting layer on each end we assume a normal metal fills the remainder of space. 
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We again assume the pair potential g(x)F(x) in each superconducting layer is a constant, 
but unlike the infinite superlattice there is no reason a priori to assume the same constant in 
each layer. We therefore assume g(x)F(x) is some constant Aj for the j th superconducting 
layer, and as before vanishes in each normal layer. This is illustrated in Fig. 3. With 
appropriate boundary conditions the Bogoliubov equations (^) are readily solved, and the 
parameters Aj are again determined by the minimization of the free energy of the system. 
This again translates into equating the spatial average of the calculated pair potential in 
the j th layer to the assumed constant Aj. These conditions now are a coupled set of N + 1 
equations, which in the limit Aj — > and AQd ^ 1 are given by 

'T c (a,b,N)~ 



ln 



7T 



16r 



-~ {K{r) - 2jyK(2jyT) + (1 + 2jy)K[(l + 2jy)r}} 



1 N A 

+ 7 E 7 ^{2{l + 3)ym + 3)yr] + 2{l-3)ym- 3 )yT} 
4 i= J+ i A i 

+ + j)y]K[l + j)yr] - [1 + (/ - j)y]K[l + (I - j)yr} 
-[-1 + (/ + j)y)K[-l + (/ + j)yr] - [-1 + j)y)K[-l + (I - j)yr]} , (21) 

where j = 0, 1, 2, . . . , N, y = 1 + b/a, and r and K(x) appear in Eq . ([T8D . As a check, if we 
consider the single layer limit with N = 0, then the system of equations (|2TD reduces to a 
single equation, which agrees with the gap equation derived before. 

We first consider the j = N th equation of the set (pT]), for which the term involving 
the summation over Ai/Aj does not contribute. This equation will determine the critical 
temperature T c (a, b, N) as a function of the thicknesses a and b of the films and of the 
number of superconducting layers present, 2N+1. In Figs. 4 and 5 we present the numerical 
solution to this equation for different values of b and N, along with the results for the infinite 
superlattice of Eq. (|19|) and for the single layer of Eq.([Tj]). Although for all values of b and 
iV the correct bulk limit as a — > oo is reached, quantitative differences with the infinite 
superlattice are obtained for smaller values of a/a c for the finite lattice. In particular, as 
b/a c increases in the finite lattice the curves do not approach that of the single layer case, 



as might be expected and is found for the infinite superlattice. 

Before concluding anything from these results, however, let us examine the other equa- 
tions of the set (|21|). For j = 0, 1, 2, . . . , iV — 1 these equations involve the term containing 
the summation over Aj/Aj, and can be used to express, for example, Ai, A 2 , . . . , Ajv in 
terms of the value of the center layer, A . In Figs. 6-9 we present the numerical solution 
to these equations at zero temperature for some representative values of b/a c and N. We 
would expect that as one moved out from the central region the superconductivity would 
become weaker in some sense, which would translate into decreasing values of the ratio 
Aj/Ao as j increases. This expectation is borne out for smaller values of b/a c for a large 
range of N, but not for larger values of b/a c . At finite temperature, illustrated in Figs. 10 
and 11, the expected behaviour is again found for smaller values of b/a c but not completely 
for the larger values, although the situation improves as the temperature increases. What 
one might conclude from this analysis is that the method used in this context seems reliable 
for thin normal films over a wide range of thicknesses of the superconducting layers, but for 
thicker normal films the method seems reasonable only for thicker superconducting layers 
(i. e., larger critical temperatures). This could indicate that the approximation of a constant 
pair potential in the superconducting layer appears reasonable in this range of validity, for 
which little spatial variation in the true potential might be expected. 

With the preceding discussion in mind we now return to a comparison of the finite 
superlattice results of the critical temperature in Figs. 4 and 5 with that of the infinite 
superlattice. We restrict this comparison to a potentially large but still finite number of 
layers in the finite system, as in the limit N — > oo our initial assumption of the independence 
of the gap parameters Aj in each layer will no longer be valid. Given this, for thicker 
superconducting layers (a/a c > 10), where the approximation of a constant pair potential 
seems reasonable over a wide range of thicknesses b/a c of the normal layer, we could conclude 
for these methods that a finite superlattice approximates well an infinite superlattice for a 
reasonably small number of layers (N ~ 10). On the other hand, for thinner superconducting 
layers (a/a c < 10), the approximation of a constant pair potential seems reasonable only for 
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thin normal layers (b/a c < 0.1). In this range we could conclude that there are noticeable 
differences between the finite lattice and the infinite superlattice in this approach, as seen 
from Fig. 4, even with a relatively large number of layers for the finite system. Qualitatively 
what is found is that, for a given thickness a/a c , the corresponding critical temperature in 
the finite lattice is smaller compared to the case of the infinite superlattice. This feature 
could be understood in terms of the form of the pair potential suggested by Fig. 3 and found 
in Figs. 6, 7, 10, and 11, where the gap parameters decrease relative to the central value as 
one moves towards the edges. If one assumes this gap parameter ratio measures in a way 
the relative strength of the binding of the Cooper pairs in each region, then from this one 
might infer that the superconductivity is in a sense "weaker" towards the edges for the finite 
potential of Fig. 3 compared to the case of the same constant gap parameter assumed over 
each superconducting region in an infinite superlattice, as in Fig. 2. Thermal excitations 
would then more readily break superconductivity in these outer regions, leading to a lower 
critical temperature for a given thickness a/a c in the finite lattice compared to the infinite 
superlattice. Such a tendency of weaker superconductivity in the finite lattice continues to 
hold except as the thickness of the superconducting layer decreases further (a/a c < 0.1); in 
the extreme limit this is a reflection of the fact that the approximation AQ<i = 2kpaVt ( i ^> 1 
used in deriving Eq. (|21~l) is breaking down. 



V. THE SINGLE LAYER REVISITED 

The results of the last section raised questions concerning the appropriateness of the 
assumption of a constant pair amplitude in a superconducting layer. In this section we will 
attempt to improve on this assumption for the case of the single layer described in Sec. II. 
Specifically, we imagine for the situation of 2N + 1 layers in Fig. 2 letting the width b of the 
normal layer approach zero. In effect we have then artificially divided a superconducting 
single layer of width 2A, as in Fig. 12, into 2N + 1 layers, each of width 2a. As in the last 
section we next assume the gap parameter in each layer is some constant Aj, and then set 
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up and solve the system of self-consistent gap equations found by a minimization of the free 
energy This again leads to a coupled set of equations, and as in the previous section the 
j = N th equation for the layer on the extreme edge will determine the transition temperature 
as a function of the film thickness a and the number of layers N. One obtains 

In (—} = - [2iVln(2iV) - (1 + 2JV) ln(l + 2N)] + - {-an(l + 2N) 
\ Qjq J 2 4 

a cos(a) - sin(a) 2N cos(2iVa) - sin(2iVa) 



a a 
+ (1 + 2N)a cos[(l + 2N)a] - sin[(l + 2N)a] 

a 

+ [2Ci(«) + a Si(a)] - (2N) [2Ci(2iVa) + 2NaSi(2Na)] 
+ (1 + 2N) [2Ci[(l + 2iV)a] + (1 + 2iV)a Si[(l + 2N)a)}} 

+ — / dxL(x) cosh" 2 (x/2r), (22) 
2r jo 



where r appears in Eq.(|18D, 

2 a T d 

a = 2Att d = — , 

na c T c 

L( x ) = 7 e + In(a;) - - [2A^ln(2A^) - (1 + 2N) ln(l + 2N)} 

~ {G{x) - 2NG{2Nx) + (1 + 2N)G[{\ + 2A^)a;]} , (23) 

G{x) is given in Eq. (p!3|) , and the Cosine and Sine integrals Ci(x) and Si(x) appear in 
Eq. (|14D . In Eq.(^) we then take a — >• and — > oo in such a way as to keep the product 
2 A = (2N + l)2a fixed. We first consider the zero temperature case, by which this limiting 
procedure will yield the critical thickness A c in this model below which superconductivity 
cannot be maintained at any temperature. We find 

1 



A c 2 T d 7e _ 1 

— = ~7^ e = ex P 

a c it l c 



(24) 



_gN(0) 

As for BCS superconductors gN(0) ranges from about 0.2 to 0.4, this critical thickness A c 
is larger by about a factor of 4 to 50 compared to the critical thickness a c obtained by the 
method described in Sec. II, where the pair potential was assumed to be constant over the 
entire superconducting layer. As well, Eq . (|24]) reinforces the trend that the critical thickness 
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would be larger for weaker BCS superconductors with smaller values of gN(0), as one might 
expect. 

As with the case of thin normal and superconducting layers of the finite lattice in the last 
section, these features of the critical thickness can be interpreted in terms of the form of the 
pair potential suggested by Fig. 12, where the gap parameters decrease relative to the central 
value as one moves towards the edges. From this one might infer that the superconductivity 
is in a sense "weaker" towards the edges for such a potential compared to the case of a 
constant gap parameter assumed over the entire superconducting region, as in Fig. 1. If 
one pictures the breaking of superconductivity as the "leaking" of Cooper pairs into the 
normal region, then one would expect that the form of the more "realistic" pair potential in 
Fig. 14 would lead to an earlier destruction of superconductivity compared to the constant 
pair potential of Fig. 1. 

One can also analyze the gap equation ( p2|) at finite temperature. For example, in the 
vicinity of zero temperature one finds 



However, at higher temperatures this method of analyzing the single layer case fails, and in 
particular gives the wrong bulk limit as A — > oo. This would indicate the assumption in this 
model of a constant pair amplitude for thick initial superconducting film layers is invalid. 



We have considered the question of when a finite sized system composed of alternating 
layers of superconducting and normal metals can effectively be thought of as a true infinite 
superlattice. Concentrating on the calculation of the critical temperature, we extended some 
methods used previously for a single superconducting layer and an infinite superlattice to 
such finite-sized systems. By examining the form of the pair potential found we argued that 
the approximation of a constant Cooper pair potential in each superconducting region seems 



Tc(A) 




(25) 



T 

- 1 - c 



VI. CONCLUSIONS 
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reasonable for thick superconducting films (a/a c > 10) over a wide range of thicknesses of 
the normal films or else for thin superconducting (a/a c < 10) and thin normal (b/a c < 0.1) 
films. The results suggest that for thick superconducting films it takes a relatively small 
number of layers (N ~ 10) using these methods to effectively make an infinite superlattice. 
However, for thin superconducting and normal films noticeable differences were found in this 
approach between the finite and infinite superlattices, even for a relatively large number of 
layers for the finite lattice. The corresponding decrease in the critical temperature of the 
finite lattice at a fixed thickness of the superconducting layer was attributed to the "weaker" 
superconductivity in the outer regions compared to the case of the constant pair potential 
assumption of the infinite superlattice. One might argue that such "weaker" zones in finite 
but large systems would be true for more realistic forms of the pair amplitude incorporating 
some spatial variation. We also applied these methods in the case of a single superconducting 
layer to construct a model which attempts to incorporate a more realistic form of the pair 
potential than the assumed constant. The model suggested that the zero temperature critical 
thickness will be enhanced with such a pair potential compared to that found by assuming a 
constant pair potential, and also exhibited the trend that weaker BCS superconductors will 
have a larger value of this thickness. These features were again attributed to the "weaker" 
superconductivity near the edges of a more realistic pair potential. 
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Figure Captions 



Fig. 1: The assumed form of the pair potential for a single superconducting layer. 

Fig. 2: The assumed form of the pair potential for an infinite superlattice. 

Fig. 3: The assumed form of the pair potential for a finite "superlattice" of 2N + 1 layers. 

Fig. 4: The critical temperature as a function of a/a c for a single layer, an infinite superlattice 
of b/a c = 0.05, and a finite "superlattice" with b/a c = 0.05 and iV = 10, 40 and 100. 

Fig. 5: The critical temperature as a function of a/a c for a single layer, an infinite superlattice 
of b/a c = 5.0, and a finite "superlattice" with b/a c = 5.0 and iV = 10, 40 and 100. 

Fig. 6: The gap parameter ratio Aj/A for the finite "superlattice" as a function of the j th 
layer at zero temperature, with N = 40 and b/a c = 0.01. 

Fig. 7: The gap parameter ratio Aj/A for the finite "superlattice" as a function of the j th 
layer at zero temperature, with N — 10 and b/a c = 0.01. 

Fig. 8: The gap parameter ratio Aj/Ao for the finite "superlattice" as a function of the j th 
layer at zero temperature, with N = 40 and b/a c = 0.5. 

Fig. 9: The gap parameter ratio Aj/A for the finite "superlattice" as a function of the j th 
layer at zero temperature, with N = 10 and b/a c = 0.5. 

Fig. 10: The gap parameter ratio Aj/A for the finite "superlattice" as a function of the j th 
layer at finite temperature T c (a, b, N)/T c = 0.283, with iV = 40 and b/a c = 0.5. 

Fig. 11: The gap parameter ratio Aj/A for the finite "superlattice" as a function of the j th 
layer at finite temperature T c (a, b, N)/T c = 0.282, with iV = 10 and b/a c = 0.5. 

Fig. 12: The assumed form of the pair potential for the single layer superconducting film of 
width 2A artificially divided into 2N + 1 layers. 
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Figure 4: b/ac=0.05 
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Figure 5: b/ac=5.0 
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Figure 6: N=40, b/ac=0.01 
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Figure 7: N=10, b/ac=0.01 
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Figure 8: N=40, b/ac=0.5 
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Figure 9: N=10, b/ac=0.5 
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Figure 10: N=40, b/ac=0.5 
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Figure 11: N=10, b/ac=0.5 
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